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Abstract. - In a cigar shaped Bose-Einstein condensate, explicit solutions of the coupled mean- 
field equations, describing defect-grey soliton dynamics are obtained, demonstrating the coexis- 
tence of grey soliton and a localized defect. Unlike the case of dark soliton, where the defect 
trapping center has vanishing superfluid density, the moving grey soliton necessarily possesses a 
finite superfluid component at the defect location. The wave vector of the impurity is controlled 
by the velocity of the grey soliton, which has an upper bound. It is found that the presence of 
the impurity lowers the speed of the grey soliton, as compared to the defect free case, where it 
can reach the sound velocity. The grey soliton's energy gets substantially modified through its 
interaction with the defect, opening up the possibility of its control through defect dynamics. 
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Quite some time back, Gross formulated the mean field 
pquations describing the macroscopic dynamics of defects 
and superfluid matter Akin to the Gross-Pitaevskii 

(GP) equation for Bose-Einstein condensate (BEC), these 
mean-field equations are well suited for describing the 
defect-BEC dynamics in a trap. On the theoretical side, 
a number of works explored the phenomenological fea- 
tures of BEC-impurity complex [3HI]- The dynamics of 
the BEC, trapped in an optical lattice, in presence of a 
localized impurity has been investigated [5]. The interac- 
tions of defect atoms in an optical lattice with a uniform 
BEC have also been studied for the understanding the de- 
phasing effect on this system j5J. Recently, Roberts and 
Rica [ID], investigated the behavior of the impurity field 
in a BEC and identified the parameter domains for the 
formation of a crystal of impurity fields and supersolid 
phases. The fact that the geometry of this complex can 
be effectively manipulated, has motivated the study of the 
BEC-defcct dynamics in lower dimensions. Cigar shaped 
BEC has received special attention, because of the pos- 
sibility of identifying exact solutions. It was found that, 
even approximate solutions in one dimension, well repre- 
sent this coupled system [TTj . 

Experimental realization of dark and bright soliton 
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P~2t417j . and soliton trains [18] has given impetus to the 
study of the interaction between defect-soliton system. 
In 1997, Konotop et al. [T5], studied the dark soliton- 
impurity complex by means of a modified adiabatic ap- 
proximation and observed differences between dark and 
bright soliton dynamics perturbed by a point defect. In 
Rcf. J3DJ , Frantzeskakis et al. , considered a static impurity 
and investigated the interaction of dark soliton with lo- 
calized impurities. It was found that the dark soliton can 
get reflected or transmitted by a repulsive impurity. Self- 
trapping of impurities was studied, for the case of both 
repulsive and attractive interactions [11] , The soliton- 
defect dynamics was also investigated by Goodman et al. 
|21| . In a recent work, it has been shown that a dissi- 
pation source can be used to generate dark, bright, gap 
and ring dark solitons, by controlling their phase and 
amplitude [22]. Very recently, Dries et al. [23], investi- 
gated the effects of impurities on collective dipole motion 
of the BEC and characterized the breakdown of superflu- 
idity of the trapped cloud in both 3D Thomas-Fermi and 
quasi-lD weakly interacting regimes. Apart from dark and 
bright solitons, GP equation supports a complex envelope 
soliton, analogous to the Bloch soliton in magnetic sys- 
tems [2~M2l)] . This grey soliton has recently been observed 
through a controlled density engineering method [2"rjll2"7] . 
Its velocity can take values from zero to the sound velocity 
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and has a dispersion very different from that of dark soli- 
ton [331[3D] . It is then natural to inquire the nature of the 
interaction between grey soliton and a point like defect. 

In a recent experiment, it was observed that, when the 
velocity of the impurity atom was made to decrease below 
the condensate sound velocity, the collisional cross section 
decreased abruptly [3D]. This indicates the better appli- 
cability of mean- field dynamics in this regime. Keeping 
this in mind, wc study here the BEC-dcfcct complex in 
the mean-field approximation, in a quasi-one dimensional 
scenario, where the defect atom moves in the condensate 
with constant velocity. As will be seen below, for this sys- 
tem both defect and soliton velocities are bounded from 
above. The general case is investigated, where defect and 
the condensate atoms have different masses. The density 
of the impurity atom is assumed to be small, so that the 
effect of the impurity on BEC excitation spectrum remains 
negligible. We find exact solutions of the coupled mean- 
field equations, describing the interaction of a grey soliton 
with a localized defect. The local minimum of the soli- 
ton in the defect location has a non-vanishing superfluid 
density. The presence of the defect lowers the grey soli- 
ton's velocity, as compared to the defect free case, where 
the maximum velocity is the sound velocity. Physically, 
the presence of the defect atom acts like a drag on the 
grey soliton, which results in a decrease of its speed from 
the sound velocity. The stability of the above solution 
is investigated using the well known Vakhitov-Kolokolov 
(VK) criteria [3T1I35] and found that the obtained solu- 
tions are stable. It is also found that the defect affects 
the energy of the grey soliton considerably, opening up 
the possibility of controlling the grey soliton through de- 
fect dynamics. We then compared these results with pure 
dark and bright soliton-defect complex, obtained from the 
mean-field equations. The mean-field solutions require the 
masses of the impurity and that of the condensate atoms 
to be same, implying that the impurity is the same atom 
in a different hypcrfinc level. Unlike the previous case, no 
restriction is found on the velocity of the soliton or that 
of the defect. 

In the model of Gross, the interaction between impurity 
and the condensate is treated in the Hartree approxima- 
tion. The defect-condensate dynamics is described by the 
following coupled equations, where the wavef unctions of 
the condensate atoms and the defect are ip a and ipb, re- 
spectively: 



where, m r 



is the reduced mass. The mass of 
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and ih 



dt 

d^ b 
dt 
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2rrib 



+ K\^ b \ 2l fa- ^a, (1) 
V 2 ^b + V ext ^ b + ^ a \ 2 ^b(2) 



Here, q = 47rh a is the strength of the atom-atom interac- 
tion and the strength of the interaction between atoms in 
the condensate and the impurity is given by, n = 2vh ° 12 , 



the condensate atoms and that of the defect are m a and 
nib, respectively, with \i being the chemical potential. The 
condensate wavefunction is normalized to the number of 
atoms N and the defect wavefunction to unity. 

For cigar shaped BEC, the three dimensional coupled 
GP equation can be transformed into quasi-one dimen- 
sion: \&i(r, t) = ipi{x 1 t)fi{y 1 z), (i = a, b). The condensate 
atoms and the impurity can be made to experience the 
same trapping potential V ex t = ^m a Lu 2 (y 2 + z 2 ), by tun- 
ing the axial frequencies of the atoms u> a and that of the 
impurity Ub, which differs by a factor of \Jm a jmb |33j . 
As is known, two different masses can bring in transverse 
separation due to gravity [33J[3I]. We assume that, this 
separation is negligible. In the reduction to quasi-one di- 
mension, taking into account the tight harmonic trap, one 
needs two different widths for the transverse Gaussian pro- 
files 
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and fb(y,z) = 



(3) 
(4) 



where a± 



and b± = 



The three dimen- 



sional coupled equations can be mapped into one dimen- 
sion by minimizing the action functional after integrating 
over the transverse degrees of freedom. In the weak cou- 
pling scenario, the coupled equations in one dimension are 
given by [551I3D] . 
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where g 



i) and k = 



(5) 
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k. We now drop 



the tilde for notational convenience and study the gen- 
eral case, where the mass of the condensate atoms and 
that of the impurity arc different. The following envelope 
profiles lead to the description of grey soliton impurity 
complex: %p a = yJ~oZe %x and ipb = ^/(Jbe' lkx ~' lu)t . It is worth 
emphasizing that for the grey soliton to exist, the defect 
should necessarily possess the plain wave component e lkx . 
Current conservation and consistency conditions yield: 



Co 



X = u(l ) and hk = rribU. 



(7) 



The soliton satisfies the boundary condition that, 
asymptotically as a — > ao, the phase variation vanishes. 
Here, oq = ^, is the equilibrium density of the conden- 
sate, u is the velocity of the soliton. The real part of Eq. 
([5]) can then be cast in a convenient form, in terms of the 
densities: 
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Fig. 1: The figure shows the variation of the grey soliton pro- 
file for different values of 8. The solid (black), dashed (blue) 
and dotted (purple) lines correspond to 8 = 0, n/8 and 7r/4, 
respectively. Inset depicts the corresponding density profiles 
of the impurity atom. The densities (|?/>i| 2 ,i = a, b) are mea- 
sured in the unit of ao. The parameter values are: aoa = 0.1, 
oj a = 2tt x 120 Hz, oj b = 2n x 122 Hz, m a = 1.44312 x 10~ 25 
Kg, m b = 1.40995 x 10 -25 Kg, a = 213a , a 12 = 99a 
(ao = 5.3 x 10 -11 m, is the Bohr radius). 
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Fig. 2: The figure shows the variation of the density of the 
impurity (Iv^fcl 2 = ff>) for different values of atom- impurity 
coupling k. The ratio m b /m a is kept fixed at 0.997 for all 
the values of g/n. The dotted (pink), dash-dotted (blue), 
dashed (purple) and solid (black) lines correspond to aia = 
50ao, 100ao, 150ao and 200ao, respectively. One can see that 
as di2, i.e., k increases, the density of the impurity increases. 
The parameter values used for this figure are same as those of 
the Fig. CO). 



+ -m a u 2 al = 0(8) 



The following densities exactly solve Eq. 
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where, b 2 = w ith cr = 

hu = p^- cos 2 



, (9) 

(10) 
and 



2£ cos 9(g / n—mt /m a ) 

kctq. The impurity atom rests 



on the minimum of the grey soliton. Notice that for ao 
to be a positive definite, — ^-J > 0. Repulsive atom- 
atom interaction leads to n < 2HI2. provided the atom- 
impurity interaction is also repulsive. If either of these 
two interactions is attractive, the localized solutions cease 
to exist. The presence of localized solitons crucially de- 
pends on the balancing between nonlinearity and disper- 
sion effects. The fact that the condensate has the self- 
interaction, apart from its interaction with the impurity, 
whereas, the impurity is devoid of self-interaction, leads 
to the dominant role of the condensate profile over the de- 
fect. It can be shown that for a significant increase in the 
value of g/n, corresponding to very weak impurity cou- 
pling, the density of the impurity is very small and has 
a marginal effect on the grey soliton. From Eq. (fTDl) . 
the following two cases (a) u = 0, i.e., static solution and 
(6) «^0, needs to be treated separately. When u = 0, 
the limit k — > leads to vanishing of the impurity wave- 
function and the grey soliton tends towards a dark soliton 
ipa — \/&o tanh(^ ). However, this limiting case leads to a 



divergence problem in the healing length ( = 



and 



y/m b Ka ' 

hence, is not admissible. For the second case, it is clear 



from Eq. (fit)]) , that for u^0, the limit k — > is unphysi- 

cal, since the density of the impurity becomes imaginary. 

2 2 

In this case, the solutions exist only when n > " " . It is 



worth mentioning that it has been recently found that, if 
the mass of the impurity atom is too heavy and strongly 
interacting, the impurity can break the condensate into 
two parts, where the mean- field theory breaks down [40] . 

Interestingly, the sound velocity in this case is found to 
be: c w = £(Jog/h. The velocity angle, also known as the 
Mach angle, is given by, 



= sin 



(11) 



where the maximum velocity of the grey soliton is 
u s = c w ^^ for 8 = 7r/2. In the present case, it is 
less than the sound velocity of the condensate. In the 
limit k — > gm a /m b , though the soliton maximum veloc- 
ity m s — > Cw, the equilibrium density of the condensate 
ao diverges. Therefore, the limit k — > gm a /m b is un- 

physical. The obtained solutions exist within the limit 

2 2 

mb(7 a < K < gm a /n%h. Superfluid density vanishes at 
9 = 0, which corresponds to a static soliton. For a mov- 
ing soliton, the superfluid density at the defect location is 
finite. 

In the presence of impurity, both the healing length and 
the maximum soliton velocity get modified. The existence 
of the localized soliton depends on the balance between 
nonlinearity and dispersion. The atom-impurity nonlin- 
ear interaction contributes to the balancing effect. As can 
be seen explicitly, the maximum velocity depends on the 
density of the condensate, interaction between the atoms 
in the condensate, as well as the interaction between the 
atoms and that of the impurity. The presence of impu- 
rity atom lowers the grey soliton's velocity. Interestingly, 
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Fig. 3: Dispersion relation of the grey soliton (solid line), in 
presence of an impurity, with the dashed line above, depicting 
the same for pure grey soliton [21]. The energy and momentum 
are measured in the units of Eo and Po, respectively. The 
parameter values, used for this figure are same as those of the 
Fig. Q. 



the effect of the impurity leads to the appearance of the 
dimcnsionlcss ratio in the maximum velocity expres- 
sion. 

The wavefunctions ip a and 4'b can also be deduced: 



ip a (x,t) = i \Z^u~ sin 6 -\- \fo§ cos 6 tanh | — - — (x — ut) 



ipb{x,t) = ibsech 



It is worth pointing out that the coupled complex en- 
velope and bright soliton solutions obtained here are of 
similar type, found earlier in the context of two compo- 
nent BEC [41] and in boson- fermion mixtures [42]. Fig. 
dU) depicts the grey soliton profiles, for various values of 
the Mach angle. Corresponding defect profiles are shown 
in the inset. It is worth observing that the localized defect 
resides at the minimum of the grey soliton, similar to the 
experimental observation in a two component BEC, when 
one component has less number of atoms [34] . Notice that, 
as 9 increases, the amplitudes of the grey soliton, as well 
as of the defect decrease. 

Fig. © shows the variation of the impurity density for 
different values of a\2- The dotted (pink), dash-dotted 
(blue), dashed (purple) and solid (black) lines correspond 
to ai2 = 50ao, 100a , 150a and 200ao, respectively. One 
can see that as a±2, i.e., k increases, the density of the 
impurity increases. For small values of k, the impurity 
density is very small, which implies that the effect of the 
impurity on grey soliton is marginal. 

The stability of the grey soliton is investigated using the 
well known criteria of Vakhitov and Kolokolov [5TJE2] ^ or 
non-linear Schrodinger type equation. It is known from 
this criterion that the solution is stable and unstable if 
H > and || < 0, respectively. When = 0, the 



solutions are found to be marginally stable. We obtained 
the exact expression of the number of atoms from the nor- 
malization condition of the condensate wavefunction, 



N = T tfjdx = — /*W (14) 
J-oo m bK V 9 

In order to obtain the stability condition, we calculate, 



dN 2h ( rribKH — u 2 rn^g 
9 V 9 



(15) 



The restriction on n gives > 0. Therefore, the obtained 
solution is found to be stable within the allowed range of 

K. 

The energy and momentum of the solitary wave and 
defect can be obtained by removing the background con- 
tribution, which ensures the convergence of the integrals, 



E, 



2m a J dx dx 



dx + ^g J (a a ~ aofdx 

+ K (a a — a )<r b dx(16) 



Using the density profiles, obtained in Eq. (TT2"j) and ((TBI . 
Eq. ([16]) yields, 



(12) 
(13) 



E a = cos 3 9 + Igafc cos 3 9-^na cos 2 9. (17) 

The three terms in the energy expression respectively 
represent, the kinetic, self-interaction and the interaction 
energy between the soliton and defect. As is evident, the 
energy of the grey soliton gets substantially modified be- 
cause of the presence of the impurity implying the possi- 
bility of controlling soliton dynamics through the defect. 
The momentum of the grey soliton remains unchanged: 
Pa = ( 7r ]^|' ^ 26* — sin 29^j . The dispersion relation is 
shown in Fig. ([BJ, which also shows the same for that of 
the pure grey soliton for comparison. The solid line de- 
picts the dispersion of the grey soliton in presence of an 
impurity, whereas, the dotted line is for pure grey soli- 
ton. Both the energy and momentum are measured in 
the unit of Eq = 7iu; a ^/o'o acr oa_L and Po = |29| . It 
is observed that for long wavelength excitations, the Lieb 
mode associated with the BEC-defect complex docs not 
differ significantly with that of the Lieb mode in defect 
free case at low momenta, and hence with the Bogoliubov 
mode. Hence, employing a single defect, one can control 
the velocity of the grey soliton, which is found to be less 
than that of the pure grey soliton. 

For the sake of completeness, we have computed the 
energy and momentum of the impurity. The energy of the 
impurity with respect to the background ctq 



E h 



h 2 P 

2nih 



h cos 2 ( 
6m fc C 2 



- k<tq cos 



(18) 
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with the corresponding momentum: Pj, = hk. 

The solutions of the coupled mean-field equations also 
admit dark and bright solitons, when both masses are 
equal (m a = m& = m): 

(19) 

and tp b (x,t) = b scch[i(a;- ut)]e lkx ~ iuj \ (20) 

with mu = hk and the healing length C = , h One 
notes that unlike the previous case, a kinematic phase 
is associated with both atoms and the impurity. The 

amplitude of the impurity is given by, b = — i— with 

V 2 C 

a o = 2g(i-i) • * s evident, for the solutions to exist, the 
atom-atom coupling must be greater than atom-impurity 
coupling. Furthermore, the solutions exist only when both 
the interactions are repulsive or attractive. As compared 
to the previous case, the velocity of the dark soliton u can 
take any finite value. 

In summary, we have obtained exact solutions for the 
grey soliton impurity complex in a cigar-shaped geometry. 
The impurity resides in the local minimum of the grey soli- 
ton, where the supcrfluid density is finite. The velocity of 
the grey soliton gets restricted, which in turn, controls the 
high momentum component associated with the impurity. 
We observed that the presence of the impurity leads to a 
reduction of grey soliton velocity, which for the defect free 
condensate can reach sound velocity. It is found that the 
impurity modifies the energy of the solitary wave, open- 
ing up the possibility to control the grey soliton dynamics 
through impurity. The obtained solutions are found to 
be stable, within the allowed range of k, as per the VK 
criteria. Pure dark or bright soliton, having a vanishing 
supcrfluid component at the defect location, are identified, 
when defect and condensate atoms have identical masses. 
Unlike the previous case, there is no restriction on the 
velocity of the grey soliton. We hope that grey soliton 
defect complex and its dynamics can be observed with 
present laboratory setup. The response of this complex 
to trap and scattering length variations [43H45] is worth 
investigating, as well as the study of this complex in an 
optical lattice [3rJHT] . 
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